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We present an investigation of the decays B" -^ A^p and B -^ ^tp^^ based on 383 x 10® 
T{4S) -^ BB decays recorded with the BABAR detector. We measure the branching fractions of 
these decays; their ratio is B{B' -> yl+p7r")/B(B'' -* A+p) = 15.4 ± 1.8 ± 0.3. The B' -^ A+p-n-- 
process exhibits an enhancement at the A^p threshold and is a laboratory for searches for excited 
charm baryon states. We observe the resonant decays B^ -^ X'c(2455)''p and B^ -^ X'c(2800)''p but 
see no evidence for B~ -^ Ec{2520)^p. This is the first observation of the decay B~ -^ I7c(2800)'^p; 
however, the mass of the observed excited X'J? state is (2846 ± 8 ± 10) MeV/c^, which is somewhat 
inconsistent with previous measurements. Finally, we examine the angular distribution of the B^ ^ 
Sci^'iSSyp decays and measure the spin of the X'c(2455)'' baryon to be 1/2, as predicted by the 
quark model. 

PACS numbers: 13.25.Hw, 13.60.Rj, 14.20.Lq 



INTRODUCTION 

Baryonic decays of B mesons, which contain a heavy 
bottom quark and a light up or down quark, provide a 
laboratory for a range of particle physics investigations: 
trends in decay rates and baryon production mechanisms; 
searches for exotic states such as pentaquarks and glue- 
balls [1, 2]; searches for excited baryon resonances; ex- 
amination of the angular distributions of B-meson de- 
cay products to determine baryon spins; and measure- 
ments of radiative baryonic B decays that could be sen- 
sitive to new physics through flavor-changing neutral cur- 
rents [3, 4]. The latter measurements rely on improving 
our theoretical understanding of baryonic B decays in 
general [5, 6]. 

The inclusive branching fraction for baryonic B decays 
is (6.8 ±0.6)% [7], and many exclusive baryonic B decay 
modes have been observed [8] . If we order the measured 
decays by Q-value: 



/ 



(1) 



where mf is the mass of each daughter in the final state 
of the B decay, we find that for each type of baryonic 
B decay, the branching fractions decrease as the Q- value 
increases. The smallest measured branching fraction is of 
the order 10^®, which also corresponds to our experimen- 
tal sensitivity for measuring these branching fractions. 
Potentially interesting i?-meson decays such as i? — > pp, 
B -^ AA, and B — > J^t^c have not yet been seen. 

Theoretical approaches to calculating baryonic B de- 
cays include pole models [9, 10], diquark models [11], 
and QCD sum rules [12, 13]. Recently, theoretical cal- 
culations have focused on pole models, where the B de- 
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cay proceeds through an intermediate 6-flavored baryon 
state, which then decays weakly into one of the final state 
baryons [14, 15]. However, it is not clear that the pole 
model is reliable for baryon poles, and the predictions 
given in the literature vary significantly. Perhaps the 
most satisfying theoretical interpretation of baryonic B 
decay rates is the qualitative one proposed by Hou and 
Soni in 2001 [16], who argue that B decays are favored if 
the baryon and antibaryon in the final-state configuration 
are close together in phase space. A consequence is that 
decay rates to two-body baryon-antibaryon final states 
are suppressed relative to rates of three-body final states 
containing the same baryon-antibaryon system plus an 
additional meson. In the three-body case, the baryon and 
antibaryon can be in the favored configuration — close to- 
gether in phase space — rather than back-to-back as in 
the two-body case. 

In this paper, we investigate the decays B'^ -^ ^tP ^^'^ 
B^ -^ yl+p7r~ [17]. We investigate baryon production in 
B decays by comparing the two-body {B'^ -^ ^tp) and 
three-body {B~ -^ yl+p7r~) decay rates directly. The 
dynamics of the baryon-antibaryon (A'^p) system in the 
three-body decay provide insight into baryon production 
mechanisms. Additionally, the B^ — > A^pn^ system is a 
laboratory for studying excited baryon states and is used 
to measure the spin of the I7c(2455)°. This is the first 
measurement of the spin of this state. 



BABAR DETECTOR AND DATA SAMPLE 

The measurements presented in this paper are based on 
383 X 10^ r(45') -^ BB decays recorded with the BaBAR 
detector [18] at the PEP-II e+e^ asymmetric-energy B 
Factory at the Stanford Linear Accelerator Center. At 
the interaction point, 9- GeV electrons collide with 3.1- 
GeV positrons at the T{AS) resonance with a ccnter-of- 
mass energy of 10.58 GeV/c^. 

Charged particle trajectories are measured by a five- 
layer silicon vertex tracker (SVT) and a 40-layer drift 
chamber (DCH) immersed in a 1.5-T axial magnetic field. 
Charged particle identification is provided by ionization 
energy (dE/dx) measurements in the SVT and DCH 
along with Cherenkov radiation detection by an inter- 



nally reflecting ring- imaging detector (DRC). 

Exclusive _B-meson decays are simulated with the 
Monte Carlo (MC) event generator EvtGen [19]. Back- 
ground continuum MC samples (e+e^ -^ qq, where 
q = u,d,s,c) are simulated using Jetset7.4 [20] to 
model generic hadronization processes. Background MC 
samples of e~^e~ -^ B^B~ and B'^B'^ are based on simu- 
lations of many exclusive B decays (also using EvtGen). 
The large samples of simulated events are generated and 
propagated through a detailed detector simulation using 
the GEANT4 simulation package [21]. 



CANDIDATE SELECTION 

We select candidates that are kinematically consistent 
with i?'^ -^ A^p and B^ -^ yl+p7r~. For the decay mode 
_B° -^ ^tPi "^6 reconstruct yl+ candidates in the pK~TT^ , 
pK^, pKgTT^TT^ , and Att^ decay modes, requiring the in- 
variant mass of each yl+ candidate to be within 10 MeV/c^ 
of the world average value [8] . For B~ -^ A'^p-K~ , we also 
reconstruct A'^ candidates in the Att'^tt^tt'^ decay mode, 
and require all of the A'^ candidates to have an invariant 
mass within 12 MeV/c^ of the world average value. 

The p, K, and tt candidates must be well-reconstructed 
in the DCH and are identified with likelihood-based par- 
ticle selectors using information from the SVT, DCH, and 
DRC. 

The Xg candidates are reconstructed from two oppo- 
sitely charged pion candidates that come from a common 
vertex; A candidates are formed by combining a proton 
candidate with an oppositely charged pion candidate that 
comes from a common vertex. The invariant mass of 
each K'^ and A candidate must be within lOMeV/c^ of 
the world average value [8] and the flight significance (de- 
fined as the flight distance from the Af vertex in the x — y 
plane divided by the measurement uncertainty) must be 
greater than 2. The mass of each K^ and A candidate is 
then constrained to the world average value [8] . 

A mass constraint is applied to all of the yl+ candi- 
dates, and all A'^ daughter tracks must come from a com- 
mon vertex. The yl+ candidates are then combined with 
an antiproton to form a i3° -^ A'^p candidate, or with 
an antiproton and a pion to form a B~ -^ yl^p7r~ candi- 
date. The daughters of each B candidate must come from 
a common vertex, and the candidate with the largest x^ 
probability in each event is selected. 

Additional background suppression is provided by in- 
formation about the topology of the events. A Fisher 
discriminant [22] is constructed based on the absolute 
value of the cosine of the angle of the B candidate mo- 
mentum vector with respect to the beam axis in the e'^e^ 
center-of-mass (CM) frame, the absolute value of the co- 
sine of the angle between the B candidate thrust axis [23] 
and the thrust axis of the rest of the event in the e+e^ 
CM frame, and the moments Lq and ^2- The quantity 



Lj is defined as J2iPi [cos^i]'', where 9i is the angle with 
respect to the B candidate thrust axis of the ith charged 
particle or neutral cluster in the rest of the event and Pi 
is its momentum. The optimal maximum value of the 
Fisher discriminant is chosen separately for each yl+ and 
B decay mode. 

Kinematic properties of _B-meson pair production at 
the T(4S') provide further background discrimination. 
We define a pair of observables, rrim and rrir, that are 
uncorrclatcd and exploit these constraints: 



2 



and 



"^^ = \l{lAtpi^-)) 



(2) 



- ms- 



The variable mm is based on the apparent recoil mass of 
the unreconstructed B meson in the event, where Qe+e- 
is the four-momentum of the e+e~ system and (Z^i+p^^-) 
is the four-momentum of the reconstructed B candidate 
after applying a mass constraint. The variable rrir is the 
difference between the unconstrained mass of the recon- 
structed B candidate and tub, the world average value 
of the mass of the B meson [8]. Signal events peak at 
TUB in mm and in rrir- This set of variables was first 
used in [24] and is chosen as an uncorrelated alternative 
to AE = Eg — ^Vs and the energy-substituted mass 



■ s — p*g (where s 



q^+ _ and the asterisk 



™es = y 4 

denotes the e+e^ rest frame), which exhibit a ^ 30% 

correlation for B^ -^ A'^pir^ . 

The event selection criteria are optimized based on 
studies of sideband data (in the region 0.10 < rrir < 
0.20 GeV/c^) and simulated signal MC samples. The data 
in a signal region (approximately ±2cr wide in m,m and 
TOr) were blinded until the selection criteria were de- 
termined and the signal extraction procedure was spec- 
ified and validated. B candidates that satisfy nim > 
5.121 GeV/c2 and [m^l < O.lOGeV/c^ are used in the 
maximum likelihood fit. 



BACKGROUNDS 

The primary source of background for i?° -^ A'^p can- 
didates is continuum e+e^ -^ qq events. Backgrounds 
due to decays such as B~ -^ AfpTT~, B° -^ A'^pn^, and 
B^ -^ S^p, S^ -^ ^t'^~ ^'"6 rejected by the criterion 
\mr\ <0.10GeV/c2. 

Approximately equal amounts of continuum e~^e^ -^ 



qq and e 



+. 



BB events make up the background 



for B^ -^ Ajpn~ events. Again, the requirement 
\mr\ < 0.10 GeV/c^ rejects most of the contributions from 
such decays as B'^ -^ yl+p7r"^7r~ and B~ ^ yl+p7r~7r''. 
Approximately 1% of the background in the fit region 
is due to these four-body events, but they do not peak 
in rUm and mr- A small peaking background is present 



from B^ -^ ^tP^ ^t ~^ ^t^^ events, especially when 
the tt" has low momentum. Based on a branching frac- 
tion measurement of the isospin partner decay B{B^ -^ 
i;c(2455)°p) = (3.7 ± 0.7 ± 0.4 ± 1.0) x 10"^ [25], where 
the uncertainties are statistical, systematic, and the un- 
certainty due to B{A^ -^ pK^n^), respectively, wc ex- 
pect 11.5 ± 2.5 peaking background events in the signal 
region for B~ — > yl+p7r~, /!+ -^ pK~Tr~^. A correction is 
applied and a systematic uncertainty is assigned to com- 
pensate for these events. 



DETECTION EFFICIENCY 

The detection efficiencies for _B° -^ A^p and B^ ^ 
A^pTi^ signal events are determined from signal MC sam- 
ples with 175, 000 to over 1, 600, 000 events in each sam- 
ple, depending on the A^ decay mode. To account for 
inaccuracies in the simulation of the detector, each MC 
event is assigned a weight based on each daughter par- 
ticle's momentum and angle. These weights are deter- 
mined from studies comparing large pure samples of pro- 
tons, kaons, and pions in MC samples and data. Small 
corrections (0.4 — 1.6%) are also appHed to account for 
tracking inefficiencies due to the displaced K^ and A ver- 
tices. These corrections depend on the K^ and A daugh- 
ter trajectories' transverse momentum and angle, and the 
distance between the beam spot and the displaced vertex. 

The detection efficiency {ei) for B^ — > A'^p signal 
events in each yl+ decay mode {I) is determined from 
the number of signal events extracted from an extended 
unbinned maximum likelihood fit to signal MC events. 
These events pass the same selection criteria as applied 
to data. The fit is performed in two dimensions, rrim 
and TTir- The probability distribution function (PDF) for 
the background consists of a threshold function [26] in 
nim multiplied by a first-order polynomial in rrir] this is 
the same as the background PDF used in the fit to the 
B^ -^ A'^p data. The signal PDF consists of a Gaussian 
in rrim multiplied by a modified asymmetric Gaussian 
with a tail parameter in rrir- The detection efficiencies in 
each /1+ decay mode are summarized in Table I. 

The detection efficiency for B~ -^ A'^p'K~ signal 
events in each yljf decay mode varies considerably across 
the Dalitz plane of the three-body decay. For reference, 
we quote the average efficiencies in Table I, but we ap- 
ply a more sophisticated treatment to these events. We 
parameterize the physical Dalitz region using the vari- 
ables cos Oh and the vI+tt" invariant mass, mA^-K- The 
helicity angle Oh is defined as the angle between the 7r~ 
and the p in the B~ rest frame. The quantity cos Oh can 
be expressed in terms of Lorentz-invariant products of 
four- vectors. We divide the kinematic region into rea- 
sonably sized bins that are uniform in cos Oh (0.2 units 
wide) and nonuniform in niA^Tr (60 — 200MeV/c^ wide). 
This choice of variables is more conducive to rectangu- 



TABLE I: Detection efficiency for B^' — > A'^p signal events, 
determined from signal Monte Carlo samples and separated 
by ylj decay mode. The numbers correspond to the efficiency 
for B" — > A^p {B^ -^ yl+pTT^), /l+ -^ fi, where /; is a given 
final state. The efficiencies quoted for the B^ -^ Atpii~ 
decays are averaged across phase space. 





Efficiency 


for A+ ^ fl 


fl 


B" - Atp 


B- 


- -> A+p-R- 


pK-n+ 


22.9% 




15.4% 


pKl 


21.6% 




14.3% 


pKln+n' 


9.6% 




5.6% 


An+ 


17.2% 




11.6% 


ATr+TT-TV+ 


- 




4.0% 



lar bins than the traditional set of Dalitz variables. The 
^TT-A^TT bins are narrower near the kinematic limits where 
the efficiency changes more rapidly and arc centered on 
expected resonances. For B~ -^ A^pn^ , yl+ -^ pK~iT^ 
near cos Oh = ^, the efficiency varies from approximately 
13% at low mA^-n, to 16% in the central niA^-n region, to 
8% at high niA^Tr- The efficiency is fairly uniform with 
respect to cos Oh, except at cos Oh ^ 1 and low mA^-n, 
where it drops to 7.4%. The other A^ decay modes ex- 
hibit similar variations in efficiency. 



SIGNAL EXTRACTION 

To extract the number of signal events in data, a two- 
dimensional {mm vs. rrir) extended unbinned maximum 
likelihood fit is performed simultaneously across /l^ de- 
cay modes. B^ -^ A^p candidates and B^ -^ A'^pn^ 
candidates are fit separately. 

The background PDF for each fit is a threshold func- 
tion [26] in rrim multiplied by a first-order polynomial in 
rrir. The shape parameter {sbkg) of the threshold func- 
tion is free but is common to all of the A'^ decay modes. 
The slope a of the first-order polynomial is allowed to 
vary independently for each yl+ decay mode. 

The signal PDF is a single Gaussian distribution in 
rUm multiplied by a single Gaussian distribution in rrir 
for B''^ -^ A'^p and multiplied by a double Gaussian dis- 
tribution in rrir for B~ —^ A'^pir" . A single Gaussian is 
sufficient to describe the signal PDF for B'^ -^ A^p be- 
cause of the small number of expected signal events. All 
of the shape parameters of the signal PDF {ssig) are free 
but are shared among the yl+ decay modes. Separate sig- 
nal (Nsig^i) and background {Nhkg,i) yields are extracted 
for each yl+ decay mode I. 

The total likelihood is the product of the likelihoods 



TABLE II: Signal yields from simultaneous fits (across A'^ 
decay modes) to B^' — > A^p and B^ -^ A^piv^ candidates. 



of the covariancc matrix V are calculated as follows: 







Ns^3 


Mode 


W -^ A+p 


B- ~* A+pTv- 


pK--K+ 


90 ± 11 


991 ± 45 


pKl 


10 ± 4 


165 ± 15 


pKl-K + TV- 


14 ± 5 


86 ± 14 


A^+ 


3 ± 3 


114 ± 13 


ATT+-K--K + 


- 


88 ± 13 



for each yl+ decay mode: 



A, 



U^iiy^ 



'l',Nsigd,Nbkg,l,Ssig,Sbkg,ai) 



(3) 



The symbol y represents the variables used in the 2-D fit, 
{■mm,mr}. 

The full simultaneous fit is validated using independent 
samples of signal MC events to simulate signal events and 
toy MC samples (generated from the background MC 
sample distribution) to represent background events in 
the fit region. For both B'^ -^ A'^p and B^ -^ A'^pn^ , 
we perform fits to 100 combined MC samples and find 
that the fit is robust and the results are unbiased. 

The results of the 2-D fits to data are shown in projec- 
tions of mm and nir for each A'^ decay mode. Figure 1 
shows the result of the fit to i?" -^ A'^p candidates and 
Fig. 2 shows the result of the fit to B~ -^ yl+p7r~ candi- 
dates. The signal yields from the fits arc summarized in 
Table II. 



BRANCHING FRACTION MEASUREMENTS 



v-/ = 



dN„dNi 



N 

E 

e=l 



fn{ye).fj{ye 



[Lk=l^kfk{ye 



(5) 



where the sum is over the N candidates. Note that in 
the calculation of the covariance matrix, the data is refit 
to the same simultaneous PDF described above, except 
that all fit parameters other than the yields are fixed to 
the values from the original fit. 

We use these sVlot weights to generate a signal or 
background distribution for any quantity that is not cor- 
related with mm or mr- The s'Plot formalism is easily 
extended to incorporate an efficiency correction for each 
candidate. Each candidate is assigned a weight of 1/e, 
where the efficiency s for an event is determined by its 
location in the cos 9h vs. m,A^TT plane. 

The branching fraction for i3° -^ A'^p for yl+ decay 
mode I is calculated as follows: 



S(B° ^ Atp)i = 



N. 



sig,l 



Nbb X £/ X 7^^ X B{A+ ^pK-n+y 



(6) 



where N^g is the number of BB events in the data sam- 
ple and TZi is the ratio of /1+ branching fraction for decay 
mode I to B{A'^ -^ pK^n^), taking care to include the 
Kg — > '!T~^TT~ and A -^ pir^ branching fractions where 
applicable. 

In order to determine the branching fraction for B^ -^ 
A'^pn^ , we take the product of the gVlot weight and effi- 
ciency weight for each candidate and sum over all of the 
candidates in the fit region. We simplify the notation by 
using sWi to denote the value of the signal sVlot weight 
for event i and include a 1% correction for the peaking 



background due to B'^ -^ ^^P, ^c 



At 



For the three-body mode B^ -^ A'^pn^ , the efficiency 
variation across the Dalitz plane requires a correction 
for each signal event in order to extract the branching 
fraction for this mode. We use the sVlot method [27] to 
calculate a weight for each event e based on the 2-D fit 
to the variables y. We have Ng = 2 species (signal and 
background) for each yl+ decay mode and define fj,k as 
the signal (j, k = 1) or background (j, k = 2) PDF. The 
sVlot weights are calculated as 



^N 



s ' nyVe 



J2j = l^njfj{ye) 
Y.kU^kIk{ye) 



(4) 



where s'Pniye) is the s'Plot weight for species n, V is 
the covariance matrix for signal and background yields, 
fj,k{ye) is the value of PDF fj^k for event e, and ye is the 
m,m and mr value for event e. The elements of the inverse 



B{B- -^ AtpTT-)i 



0.99 X (E — 



N 



BB 



X TZi X B{At ^pK-TT+y 



(7) 



The contribution from the peaking background is esti- 
mated using the yl+ -^ pK^ir^ decay mode. Since the 
overall branching fraction for the peaking background 
contribution is the same regardless of /l^ decay mode, 
it is applied as a proportional correction. The measure- 
ments for B{B'^ -^ A'^p) and B{B~ — > yl+p7r~) for each 
yl+ decay mode are summarized in Table III. 

The BLUE (Best Linear Unbiased Estimate) technique 
is used as described in Ref. [28] to combine the corre- 
lated branching fraction measurements for different /!+ 
decay modes. The purpose of the method is to obtain 
an estimate x that is a linear combination of t individual 
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FIG. 1: Projections of irim (left) and rrir (right) in data for S" — > /l^p candidates, separated by A^ decay mode: (a, 6) are 
A'^ ^ pK~n^ , (c, d) are A^ ^pKs, (e, /) are yljt -^ pKsTv'^'!v~ , and (<?, ft) are yljt — > tItt^. The TTim projections {a,c,e,g) are 
for Inirl < 0.030 GeV/c^ and the TTir projections {b,d,f,h) are for nim > 5.27GeV/c^. The solid curves correspond to the PDF 
from the simultaneous 2-D fit to candidates for the four At decay modes, and the dashed curves represent the background 
component of the PDF. 



measurements (x/), is unbiased, and has the minimum 
possible variance a^. The estimate x is defined by 



E 



aixi. 



(8) 



The condition J^i «/ = 1 ensures that the method is un- 
biased. Each coefficient a/ is a constant weight for mea- 
surement xi and is not necessarily positive. The set of 
coefficients ct (a vector with t elements) is determined by 



UtE-iU' 



(9) 



where U is a i-component vector whose elements are all 
1 (Ut is its transpose) and E is the (i x t) error matrix. 
The diagonal elements of E are the individual variances, 
af. The off-diagonal elements are the covariances be- 
tween measurements (raiai', where r is the correlation 



between measurements / and I'). The error matrices add 
linearly, so we define E = Estat + 'Esyst- Estat includes 
the uncertainties in the fit yields and the correlations 
between yields from the simultaneous fit result. Egyst in- 
cludes the systematic uncertainties that are described in 
the next section. Overall multiplicative constants {Ngj^ 
and B{At -^ pK~7T'^)) that are common to all the mea- 
surements and their uncertainties are not included in the 
BLUE method. 



The solutions for a are 
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), (10) 



^ 




100- 



f^^'i^ k '^'^x*^ ' ^^^, 



« 


80 


o 




oo 


60 


o 




o 


40 


o 




"^ 


20 


r/1 




-^^ 






80 


> 




w 


60 




40 




20 



100 



i^ ^A^4^ 



M^^*T?+M*r 



M^^"^w^ 




5.15 



5.20 



5.25 



5.30 
m„, (GeV/c^) 



300 
200 
100 



60- 



40 



20 

> 

O 40 

oo 
o 
o 




*^^^ 






20 T 



40- 



* ^4fM^HJ 








'i^ '^'^^ • A: 



-0.10 



-0.05 



0.00 



0.05 



0.10 



m^ (GeV/cO 



FIG. 2: Projections of mm. (left) and nir (right) in data for B~ — > yl^pyr" candidates, separated by A'^ decay mode: (a,fe) are 
At -^ pK~Tv~^, {c,d) are At -^ pKl, (e,/) are At —> pK^-k^-k^ , {g,h) are At -^ A-k'^ , and (i,j) are /!+ ^ yl7r+7r~7r+. The 
mm projections {a,c,e,g,i) are for \m,r\ < 0.030 GeV/c^ and the rrir projections {b,d, f,h,j) are for m™, > 5.27GeV/c^. The 
soUd curves correspond to the PDF from the simultaneous 2-D fit to candidates for the five At decay modes, and the dashed 
curves represent the background component of the PDF. 



and 

B- -^ AtpTT^ : 

aT = ( 0.913 0.043 -0.003 0.029 0.018 ) , (11) 

where cxt is the transpose of a. The order of the co- 
efficients corresponds to the order of At decay modes 
presented in Table III. 

We calculate the best estimate x according to Eqn. 8 
and divide this quantity by N^g and B{At -^ pK^n^). 



We calculate the variance of x 



a = OLT E a. 



(12) 



Since the error matrices add linearly, we quote separate 
statistical and systematic uncertainties. The statisti- 
cal and systematic uncertainties in iV^g- are added in 
quadrature with the statistical and systematic a results, 
respectively. The combined branching fraction measure- 



10 



ments arc thus 

S(S° ^ Atp) = 

(1.89 ± 0.21 ± 0.06 ± 0.49) x 10"^ 
B{B-^AtpTT-)^ 

(3.38 ± 0.12 ± 0.12 ± 0.88) X 10"^, (13) 

where the uncertainties arc statistical, systematic, and 
the uncertainty in B{A^ ^ pK~Tr~^), respectively. 

We use the same procedure to determine the branching 
ratio B{B- -^ A+pTr-)/B{B" -^ A+p): 



B{B- -^ A+piT- 
B{B" ^ A+p) 



= 15.4±1.8±0.3. 



(14) 



In the branching ratio, many of the systematic uncer- 
tainties, including the dominant B{A+ -^ pK~TT+) un- 
certainty, cancel. 



SYSTEMATIC UNCERTAINTIES IN THE 
BRANCHING FRACTIONS 

The uncertainties in the B{B'^ -^ ^tp) ^^^ B{B^ -^ 
A+pn^) measurements are dominated by the uncertainty 
in the A+ —^ pK~iT+ branching fraction, and then by 
the uncertainties in the A+ branching ratios (compared 
to A+ ^pK-'K+) [8]. 

The systematic uncertainties for each A+ decay mode 
are summarized in Tables IV and V. The systematic 
uncertainty in the number of BB events is 1.1%. 

The uncertainty in the efficiency determination is due 
to MC statistics, charged particle tracking, and particle 
identification. For _B° -^ A+p, the uncertainty due to MC 
statistics is (0.4 — 0.6)%. For B^ -^ A+p'K~ , we calculate 
the uncertainty due to MC statistics by independently 
varying the number of reconstructed signal MC events 
in each cos 0h,mA^T7 bin according to a Poisson distribu- 
tion (ensuring that the data events in the same bin are 
correlated). The resulting uncertainty is (0.6 — 3.1)%. 

The tracking efficiency systematic uncertainties are de- 
termined from two separate studies. In the first study, 
T+T~ candidates are selected, in which one r candidate 
decays to leptons and the other decays to more than one 
hadron plus a neutrino. Events are selected if one lepton 
and at least two charged hadrons are reconstructed. The 
efficiency is then measured for reconstructing the third 
charged particle for the hadronic decay. From this study 
there is a (0.38 — 0.45)% uncertainty in the tracking ef- 
ficiency per charged particle. In the second method, a 
charged particle trajectory is reconstructed in the SVT, 
and then the efficiency for finding the corresponding tra- 
jectory in the DCH is measured. For the latter study, 
the uncertainties range from 0.21% to 1.18% depending 
on the A+ decay mode. The systematic uncertainties 
determined in the two studies are added in quadrature. 



The systematic uncertainty for charged particle iden- 
tification is a measure of how well the corrections ap- 
plied to the events in the signal MC sample for the 
efficiency determination describe the i?° -^ A+p and 
B~ — > A+pTT~ decay modes. The corrections are deter- 
mined from control MC and data samples (a yl ^ pn^ 
sample for protons and D*+ -^ 0^11+, D'^ -^ K^ti+ 
samples for pions and kaons). The efficiency as a func- 
tion of momentum and angle for the B^ -^ A+p and 
_B~ -^ A+p-iT~ signal MC samples and the control MC 
samples agree to within (0.8 — 3.5)% for different ranges 
of momentum and angle. 

Systematic uncertainties due to the fit procedure are 
also considered. The dominant fit uncertainty is due 
to the threshold function parameter in the background 
PDF. The fit validation study showed that this param- 
eter must be shared among A+ decay modes to ensure 
fit robustness. We allowed this parameter to vary inde- 
pendently among A+ decay modes and repeated the fit 
to data; the difference between the fit results is taken as 
a systematic uncertainty of (1 — 24)% depending on the 
A+ decay mode. A peaking background due to misre- 
constructed -B'^ -^ ^tP-i ^t ^ ^t^'^ events is present 
at the level of 1.0%. We assign a systematic uncertainty 
of 0.5% to account for the uncertainty in B(B'^ -^ ^cP)- 
The nominal endpoint in the fit to mm is 5289.0 MeV/c^; 
we vary this by ±0.5MeV/c^, resulting in a systematic 
uncertainty of (0.2 — 1.5)%. 



A+p THRESHOLD ENHANCEMENT 

The kinematic features and resonances in _B^ -^ 
A+pT:~ are investigated through examination of the 2-D 
Dalitz plot (wp^ vs. m^^^) and its projections (the res- 
onances will be discussed in the next section) . Using the 
aVlot formalism, we project the events in the {mm, m,r} 
fit region using the signal sVlot weights and background 
sVlot weights along with the efficiency corrections. This 
method allows us to project only the features of the signal 
events, while taking the efficiency variations into account. 
Figure 3 shows the gVlot weights for m^^ vs. m,\^^. Note 
that the negative bins are suppressed in the 2-D Dalitz 
plot. 

We project the events in the fit region onto the niA^p 
axis with signal s'Plot weights and efficiency corrections 
to study the enhancement at threshold in the baryon- 
antibaryon mass distribution. This enhancement can be 
seen in B~ — > A+pTr~ decays as a peak in m,A^p near 
the kinematic threshold, m^^^^ = 3224.8 MeV/c^. We di- 
vide the normalized ruA^p distribution by the expectation 
from three-body phase-space; the resulting distribution 
is shown in Figure 4. An enhancement is clearly visi- 
ble near threshold. The observation of this enhancement 
is consistent with baryon-antibaryon threshold enhance- 
ments as seen in other decay modes such as _B -^ PpK, 
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TABLE III: Summary of the individual and combined branching fraction measurements for B^' 
The uncertainties are statistical, systematic, and due to B{A'^ -^ pK^n'^), respectively. 



A^p and B 
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At 


^pKl 




At 


^pKU+ 


■K 


At 


~.A^+ 




At 


-^ ylTT + TT" 


TT 


combined 





(2.05 ± 0.25 ± 0.05 ± 0.53) x 10" 

(1.49 ± 0.60 ± 0.17 ± 0.39) x 10" 

(4.33 ± 1.55 ± 0.57 ± 1.13) x 10" 

(0.71 ± 0.71 ± 0.18 ± 0.18) X 10" 

(1.89 ± 0.21 ± 0.06 ± 0.49) x 10" 



(3.38 ± 0.13 ± 0.11 ± 0.88) x 10"* 
(3.82 ± 0.35 ± 0.38 ± 0.99) x 10"* 
(4.58 ± 0.70 ± 0.66 ± 1.19) x 10"* 
(3.98 ± 0.45 ± 0.39 ± 1.03) x 10"* 
(3.49 ± 0.51 ± 0.38 ± 0.91) x 10"* 
(3.38 ± 0.12 ± 0.12 ± 0.88) x 10"* 
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FIG. 3; Dalitz plot of m^^ vs. m\^^. Each event is efficiency-corrected and given a signal (left) or background (right) sVlot 
weight. Note that the density scales on the left- and right-hand plots are different, and negative bins are suppressed. 



B -^ Dpp{Tr), or in contimmin e^ 



■PPl [1, 2]. 



RESONANT SUBSTRUCTURE OF B" 



A+pn- 



We also project the events in the fit region onto the 
niA^TT axis with signal gVlot weights and efficiency cor- 
rections to study resonances in vriA^-n- We perform 1-D 
binned y^ fits to discriminate between resonant (B" ^ 
£'" p) and nonresonant (i?" -^ Atpn^) signal events. 
In each binned x^ fit? the PDF is numerically integrated 
over each (variable-sized) bin and the following quantity 
is minimized: 



rib 






J nr / nr j CtTfli 



Y, 



(15) 



where Vsig is the resonant signal PDF, T'm. is the nonres- 
onant signal PDF, Fgjg is the expected yield of weighted 
resonant signal events, and Fnr is the expected yield of 
weighted nonresonant signal events. We assume the am- 
plitude and phase of the non-resonant B~ -^ Atp'!T~ 
contribution is constant over the Dalitz plot, and does 



not interfere with the resonant contributions. The range 
of the integral over the quantity drrii takes into account 
the variable bin width, Yi is the number of weighted 
data events, and cFi is the uncertainty in Yi. Variable 
bin widths are used to ensure that there are a sufficient 
number of signal events in each ttia^tt bin so that the es- 
timated uncertainty is valid. This is especially important 
in the non-resonant sideband regions. Bin widths in the 
signal regions are chosen to have sufficient granularity 
throughout the resonant peaks. 

The Z'c(2455)° and Z'c(2520)° are well-established res- 
onances that decay to tI^tt". A third Ec resonance, the 
Z'c(2800)°, was reported by the Belle Collaboration in 
2005 [29] along with its isospin partners i7c(2800)+ and 
Z'c(2800)++. These resonances were observed in con- 
tinuum (e+e" -^ cc) AcTT events. The Sd'^SOO)'^ reso- 
nance was fit with a D-wavc Brcit-Wigner distribution 
and the mass difference Am = ms^(28oo)'> ~ "^/i+ ~ 
(515 ± 3 t g) MeV/c^ was measured, which corresponds 
to an absolute mass of (2802 t 7) McV/c^ [8]. The nat- 
ural width of the resonance is (61^; ^g)MeV [29]. We 
search for evidence of all three resonances. 



12 



TABLE IV: Summary of the contributions to the relative systematic uncertainty in B{B —> Atp) for each At decay mode. 
The total for each mode is determined by adding the uncertainty from each source in quadrature. The fractional statistical 
uncertainty in the fit yield for each mode is provided for comparison. 



Source 

At 


-^pK^Tv'^ 


-+ Atp Systematic Uncertainty 
At -^ pK°s At -^ pKU+-k- yl; 


t ^AlT + 


BE events 


1.1% 


1.1% 




1.1% 




1.1% 


TZi 


- 


8.5% 




11.8% 




8.9% 


MC statistics 


0.4% 


0.6% 




0.6% 




0.4% 


Tracking 


1.7% 


1.9% 




2.8% 




1.7% 


Displaced Vertices 


- 


1.1% 




1.1% 




1.1% 


Particle Identification 


1.5% 


2.1% 




1.7% 




1.6% 


Fitting 


0.9% 


7.0% 




4.9% 




24.2% 


Total 


3% 


12% 




13% 




26% 


Statistical 


12% 


40% 




36% 




100% 



TABLE V: Summary of the contributions to the relative systematic uncertainty in B{B^ — » AtpTr^) for each At decay mode. 
The total for each mode is determined by adding the uncertainty from each source in quadrature. The fractional statistical 
uncertainty in the fit yield for each mode is provided for comparison. 



Source 

At 


— * pK 


B' -+ Atpiv' 
"7r+ At -^ pK^s At - 


Systematic Uncertainty 

^ pKl-K + TX- At -^ AtX+ At ' 


-^ ATX + -K-TX + 


BB events 


1.1% 


1.1% 


1.1% 


1.1% 


1.1% 


111 


- 


8.5% 


11.8% 


8.9% 


6.1% 


MC statistics 


0.6% 


2.1% 


3.1% 


2.0% 


3.0% 


Tracking 


2.6% 


2.3% 


3.2% 


2.1% 


2.9% 


Displaced Vertices 


- 


1.1% 


1.1% 


1.1% 


1.1% 


Particle Identification 


0.8% 


1.8% 


2.5% 


1.8% 


3.5% 


Fitting 


1.6% 


3.2% 


6.5% 


2.5% 


2.4% 


Total 


3.4% 


9.9% 


14.5% 


10.0% 


8.7% 


Statistical 


4.5% 


9.1% 


16.3% 


11.4% 


14.8% 
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FIG. 4: Projection of the amplitude squared (1^4^) vs. niAap 
for B^ -^ A^pn^ decays near threshold. Candidates are 
efficiency-corrected, weighted using the sVlot technique, and 
corrected according to three-body phase space. 



A significant i7c(2455)'^ signal is seen near threshold; 
see Fig. 3 and Fig. 5. We construct a resonant signal PDF 
from a non-relativistic S'-wave Breit-Wigner distribution 
convolved with the sum of two Gaussian distributions (to 
form a "Voigtian" distribution). This quantity is multi- 
plied by a phase-space function. The mass and (constant) 
width of the resonance in the Breit-Wigner PDF are free 
in the fit. The resolution (described by the two Gaus- 
sians) is fixed; it is determined from a B~ — > Z'c(2455)"p, 
i:c(2455)'' -^ A+TT-, A+ -^ pK-7T+ signal MC sample 
by comparing the measured Z'c(2455)° mass to the true 
i3'c(2455)'' mass for each candidate. The two-body phase- 
space function goes to zero at the kinematic threshold for 
A^tt: 2426.03 MeV/c2. The non-resonant signal PDF is a 
threshold function [26]; the threshold is set to the kine- 
matic threshold for A^tt. The shape parameter for the 
threshold function is fixed from a fit to a non-resonant 
B'^ -^ yl+p7r~ signal MC sample. The weighted data 
points are shown in Fig. 5 with the fit overlaid; the fit 
results are summarized in Table VI. The average effi- 
ciency for yl+ -^ pK^TT^ signal MC events in this region 
is 14.1%. 

No significant signal is seen in the region of the 
Z'c(2520)*'; see Fig. 6. We perform a fit using a rela- 
tivistic D-wave Breit-Wigner distribution with a mass- 
dependent width to describe the resonant signal PDF, 
fixing the resonance mass and the width to the world 
average values [8]: mn = (2518.0 ± 0.5) MeV/c^ and 
Tr = (16.1 ± 2.1) McV. The non-resonant signal PDF 
is a first-order polynomial. We obtain Ysig = 27 ± 69 
events; the fit result is shown in Fig. 6. The average effi- 
ciency for yl+ -^ pK^TT^ signal MC events in this region 
is 15.4%. 

In the niA tt distribution, we also observe an excited 
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FIG. 5; (a) Projection of mAcn showing the X'c(2455) reso- 
nance. Events are efficiency corrected and weighted using the 
sVlot technique, and the result of a binned x^ fit to a Voigtian 
signal plus a threshold function background is overlaid. The 
variable bin sizes range from 1 to 7MeV/c^. (b) The same fit 
result and data is shown on a smaller vertical scale to show 
the behavior of the PDF at threshold. 



TABLE VI: Fit resuhs for B' -^ rc(2455)°p. Kig is the 
efficiency-corrected resonant signal yield in the fit range. Sys- 
tematic uncertainties from the fit to rrim vs. rrir are not 
included in the yield. The world average values from the 
Particle Data Group (PDG) of the mass and width of the 
X'c(2455)'' are included for comparison [8]. 



Fit Parameter 



Value 



PDG Value [8] 



Y,i 



1522 ± 149 



mR ( GeV/c^) 2.4540 ± 0.0002 2.4538 ± 0.0002 
Tr (MeV) 2.6 ± 0.5 2.2 ± 0.4 
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2.48 2.50 2.52 2.54 2.56 2.58 2.60 2.62 2.64 



m^ „ (GeV/c^) 

FIG. 6: Projection of niA^TT in the region of the X'c(2520)'^ 
resonance. Events are efficiency corrected and weighted using 
the sVlot technique, and the result of a binned x^ fit to a rel- 
ativistic _D-wave Breit-Wigner signal with a mass-dependent 
width plus a linear background is overlaid. The bin size is 
5MeV/c^. No significant signal is seen. 



S^ state that is higher in mass than the Z'c(2520)*'. We 
construct a resonant signal PDF from a relativistic Breit- 
Wigner distribution with a mass-dependent width: 



r(g) = T;,. - 

\qR 



2LH 



/ mR 



\mA, 



B'L{q,qR? 



(16) 



The quantity L is the angular momentum (L = 0, 1, 2 is 
S'-wave, P-wave, D-wave, respectively), q is the momen- 
tum of the A^ (which is equal to the momentum of the 
7r~) in the excited I7j? rest frame, and qn is the value 
of q when ruA^T^ = tur. In Eqn. 16, B'^{q^qR) is the 
Blatt-Weisskopf barrier factor [8] : 

Kiq^qR) = 1' 



(17) 



where we define a constant impact parameter d = 1 f m 
(the approximate radius of a baryon) , which corresponds 
to 5.1 GeV~^. Blatt-Weisskopf barrier factors are weights 
that account for the fact that the maximum angular mo- 
mentum (L) in a strong decay is limited by the linear 
momentum (q). Since the resonance is quite wide, we do 
not need to include a resolution function in the resonant 
signal PDF. The two fit parameters {niR and Tr) of the 
S^ are free in the fit. The non- resonant signal PDF is a 
first-order polynomial. 

From the 1-D binned x^ fit, we obtain Fgjg = 1449±284 
efficiency-corrected events for the excited E^ state. We 




TABLE VII: Fit results for the excited X'J? resonance. Ksig is 
the efficiency-corrected resonant signal yield in the fit range. 
Systematic uncertainties from the fit to nim vs. nir are not 
included in the yield. The world average values from the 
Particle Data Group (PDG) of the mass and width of the 
X'c(2800)° are included for comparison [8]. 



Fit Parameter 



Value 



PDG Value fSl 



Y ■ 

niR ( GeV/c 
Tr (MeV) 



1449 ± 284 
2.846 ± 0.008 
86 



+ 33 
22 



2.802 
61 : 



- 0.004 
■ 0.007 

28 

18 




2.5 2.6 2.7 2.8 2.9 3.0 3.1 3.2 3.3 



m^^„ (GeV/c^) 

FIG. 7: Projection of niA^rr showing an excited E^^ resonance. 
Events are efficiency corrected and weighted using the aVlot 
technique. The result of a binned x^ fit to a relativistic D- 
wave Breit-Wigner signal with a mass-dependent width plus 
a linear background is overlaid. The variable bin sizes range 
from 15 to 40 MeV/ c^ 



choose L = 2 for the nominal fit, but investigate L = 0, 1 
as well. The average efficiency for yl+ — *■ pK~Tr~^ signal 
MC events in this region is 16.3%. The fit results are 
summarized in Table VII and shown in Fig. 7. The x^ 
from the fit is 37 with 31 degrees of freedom (DOF). If 
the signal yield is fixed to zero and the mean and width 
are fixed to the central values from the nominal fit, the 
resulting x^ is 78 with 34 DOF. The significance can be 
calculated from A^^ — 40.9, which is equivalent to 5.8(t 
for the joint estimation of three parameters. 

The measured width of this state (86 1 22 ± 12) MeV 
is consistent with the width of the Ec{2800)'~' measured 
by Belle [29]. However, the measured mass of this ex- 
cited i:° is (2846 ± 8 ± 10) MeV/c^, which is 40MeV/c2 
and 3ct higher (assuming Gaussian statistics) than Belle's 
measured mass for the i7c(2800)°. 

We evaluate systematic uncertainties for the Z'c(2455)° 
and the excited S^ yields, masses, and widths by modify- 
ing the binning, the resonant signal PDF shape, and the 



15 



TABLE VIII: Systematic uncertainties for Ysig, m_B (in 
MeV/c^), and r^ (in MeV) for the X'c(2455)" and excited 
X"!? resonances. 



Systematic Source 



X'c(2455)" excited El 
Vu "Ksig mu Tfl 



Y ■ 



Resonant Signal PDF - - 5.9% - ± 7 

Non-resonant Signal PDF - - 1.2% - ± 2 

Binning 6.9% ±0.3 20% ±10 ±10 



Total 



6.9% ±0.3 21% ±10 ±12 



non-resonant signal PDF shape. Changing the variable 
bin sizes leads to the dominant systematic uncertainty in 
the masses, widths, and yields of both resonances. For 
the Z'c(2455)'^, the bin width in the peak region was de- 
creased from the nominal IMeV/c^ to 0.5 MeV/ c^; for 
the excited 17^, the bin width was varied from 10 to 
20 MeV/c^ compared to the nominal 15 MeV/c^. For both 
resonances, an S'-wave and a P-wave relativistic Brcit- 
Wigner (without a resolution function) was used instead 
of the nominal resonant signal PDF. The (fixed) non- 
resonant threshold parameter for the £"0(2455)° was var- 
ied by ilcr. A second-order polynomial was used (in- 
stead of a first-order polynomial) for the excited Z'J? non- 
resonant PDF. A summary of the systematic uncertain- 
ties for Vsig, JTiij, and V^. are summarized in Table VIII. 
The significance is recalculated following each of the 
variations used to evaluate the systematic uncertain- 
ties in the excited Z"" resonance parameters. The re- 
sulting significance (including systematics) is 5.2(T. A 
cross-check is performed to make sure the Z'c(2800)" sig- 
nal is not the result of interference with a A(1232)++, 
for example (although no significant A(1232)++ signal 



is seen in the 



''PTT 



distribution). The fit is performed 



again in the Z'c(2800) mass region for candidates with 
niipT^ > 1.5GeV/c^. We obtain 1329 ± 230 resonant signal 
events (compared to 1449 ± 284 events for the nominal 
fit) and a consistent mass and width. 

An additional cross-check is performed to investi- 
gate whether there are appropriate fractions of resonant 
Z'c(2800)'' events in different yl+ decay modes. This is 
accomplished by dividing the s'P lot-weighted, efficiency- 
corrected data into two samples according to the yl+ de- 
cay mode. Note that this cross-check neglects statistical 
correlations from the combined rrir vs. mm fit (less than 
15%) among the yl+ decay modes. A binned x^ fit to only 
yl+ -^ pK"^^ candidates gives l^ig = 776 ± 160, com- 
pared to 6463 ± 241 total non-resonant B^ -^ A'^pn^ , 
A+ -> pK-TT+ events ((12 ± 3)%). A binned x^ fit to 
a combined sample of yl+ -^ pKg, ^J — > pK g7r~^Tr~ , 
yl+ -^ AiT^, and yl+ -^ Att^tt^tt^ candidates gives 
l^ig = 530 ± 177 compared to 5956 ± 431 non-resonant 
events ((9 ± 3)%). (In order for this fit to converge, mji 



and Tn were fixed to their nominal values.) The fractions 
are consistent in the two samples and the total (1306±239 
events) is consistent with the nominal fit result within 
uncertainties. 

We have also investigated the possibility that there 
are two resonances in the mass range shown in Figure 7. 
However, there is no evidence for two distinct vertical 
bands in this region in the B^ -^ yl+p7r~ Dalitz plot 
(Figure 3), and we do not obtain a statistically significant 
fit to two resonances. 

In order to measure the fraction of B^ -^ A'^pi:^ de- 
cays that proceed through intermediate Sc resonances, 
we assume that the contribution from each A'^ decay 
mode for events in the Sc regions is equal to the measured 
contribution from each A'^ decay mode in all (resonant 
and non-resonant) B^ -^ A'^pir^ events. We set a 90% 
C.L. upper limit on B^ -^ Uc{2520)^p that includes sys- 
tematic uncertainties and corresponds to 109 events. The 
measured fractions or upper limits of B^ -^ A^p'K~ de- 
cays that proceed through an intermediate Sc resonance 
are 



B{B- -^ Sc{2Abhfp) 



B{B- 
B{B- - 



i:c(2800)"p) 



B{B- 
B(B- - 



-> Ajpn-) 
i:c(2520)0p) 



B{B- -^ Atpii- 



(12.3 ±1.2 ±0.8) X 10" 



(11.7±2.3±2.4) X 10" 



<0.9x 10^2(90% C.L.). 



(18) 



Therefore approximately 1/4 of _B — > A'^pir decays 
proceed through a known intermediate Sc resonance. 



MEASUREMENT OF THE 1:^(2455)'' SPIN 

The -£'c(2455)° is the lowest mass Sc state. In the 
quark model, it is expected to have J^ = 2 ? where J is 
the spin and P is the parity. In this section, we provide 
a quantitative evaluation of the spin-1/2 and spin-3/2 
hypotheses for the £"£(2455)" baryon. 

We determine the spin of the £"£(2455)° through 
an angular analysis of the decay B^ -^ Z'c(2455)°p, 
Z'c(2455)° -^ AfTT~. We define a hehcity angle 9^ as the 
angle between the momentum vector of the A'^ and the 
momentum vector of the recoiling _B-daughter p in the 
rest frame of the Z'c(2455)°. If we assume J(yl+) = 1/2, 
the angular distributions for the spin-1/2 and spin-3/2 
hypotheses for the I7c(2455)" are 






dN 



d cos Oil 

dN 
d cos dh 



ex 1 



(19) 



ex l±3cos^0,i. 



These are the ideal distributions; the measured angular 
distributions will be somewhat degraded due to nonuni- 
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form detector efficiencies, finite experimental resolution 
for measuring Oh, and background contamination. We 
estimate the effects of inefficiencies and background con- 
tamination by performing parameterized MC studies to 
quantify the decrease in sensitivity to discriminate be- 
tween possible spin values. 

If we select from a ±2(7 signal region in mm and Wr, 
without sVlot weights or efficiency corrections, there are 
127 events in the i7c(2455)° signal region and 27 events 
in the Sc{24:55)'^ background regions (2.430 < mA^TT < 
2.445 GeV/c2 and 2.463 < niA,^ < 2.478 GeV/c^). Scal- 
ing the number of events in the background region by 
the ratio of the total width of the background regions 
compared to the width of the signal region, we expect 
7.2 ±1.3 background events in the signal region. 

We measure the finite experimental resolution of cos 6h 
by comparing the measured value of cos 6h to the true 
value of cos 9h in B^ -^ A'^pir^ , yl+ -^ pK^n^ events in 
a signal MC sample. The maximum root mean square of 
the measured value of cos Oh minus the true value of cos Oh 
in the Z'c(2455)° signal region determines the helicity 
angle resolution (j(cos Oh) < 0.03. Therefore the finite 
experimental resolution is small compared to any features 
in the spin-1/2 or spin-3/2 angular distributions. 

We investigate the discrimination power between spin 
hypotheses using parameterized MC studies. In general, 
the log hkelihood is computed as ln£ = ^jWiln(2/i), 
where yi is the probability density for observing event i. 
The weight Wi for the MC studies is Wi = Si, where Ei 
is the efficiency for event i. We compute a log likelihood 
for each hypothesis: 



ln/:(l/2)=^w;,lnj 

i 

ln£(3/2) ^^^wjln 



{l + 3 cos^ Oh., 



(20) 



The shape of the cos Oh distribution for background 
events is estimated from the shape of the helicity distri- 
bution for events in the background regions. The helicity 
distribution for these events is illustrated in Fig. 8 as a 
non-parametric PDF (a histogram). This PDF is used to 
generate the number of background events in the signal 
region with a Poisson uncertainty (7.2 ± 2.7). The total 
number of events in the sample is fixed to 127, so the 
background effectively dilutes the signal distribution. 

We then generate 500 samples (127 events each) and 
compute the likelihood £ that each generated distribu- 
tion is uniform in cos Oh (spin-1/2) or distributed as 
1 + 3 cos^ Oh (spin-3/2). We define the quantity A ln£ = 
ln£(l/2) — ln£(3/2). Figure 9 shows the distribution 
Aln>C for events generated with each hypothesis. The 
dashed histogram (negative values of A In £) corresponds 
to samples generated according to the spin-3/2 hypoth- 
esis, and the solid histogram (positive values of Aln£) 
corresponds to samples generated according to the spin- 
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FIG. 8: Helicity angle distribution for the combined sam- 
ple of background and non-resonant signal events, a non- 
parameteric PDF (line). The shaded region indicates the 
Poisson uncertainty in the distribution. 



1/2 hypothesis. For each distribution, the separation 
from zero illustrates how well we can discriminate be- 
tween hypotheses given 127 signal events. 



The helicity angle distribution for events in the signal 
region around the Z'c(2455)'^ is shown in Figure 10. The 
points are efficiency-corrected. Functions corresponding 
to the spin-1/2 (solid) and spin-3/2 (dashed) hypotheses 
are overlaid. We compute the difference in log likelihood 
between the hypotheses: Aln£ = +19.2. Wc indicate 
the value of A In/! in data with a vertical line in Figure 9. 
The observed value of A In £ is consistent with the spin- 
1/2 hypothesis and excludes the spin-3/2 hypothesis at 
the > 4(7 level. 



The ideal angular distributions for the Z'c(2455)° 
stated in Eqn. 19 are also applicable for the excited S^ 
resonance. Unlike the narrow Z'c(2455)° resonance near 
threshold, the excited E^ is much wider and therefore its 
angular distribution is extremely contaminated by the 
non-resonant signal events underneath the signal. We 
perform a non-resonant sideband subtraction to extract 
the helicity angle distribution of the excited S^, but are 
limited by the number of signal events available. An 
examination of this distribution is somewhat consistent 
with a J = 1/2 hypothesis, but no conclusive statement 
can be made about the spin of the observed excited S^. 
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A ln(L) 

FIG. 9: Distribution of A In £ = in £(1/2) - In £(3/2) for 
signal events generated with a uniform distribution in cos Oh 
(solid histogram, positive values) and a 1 + 3 cos^ 6h distribu- 
tion (dashed histogram, negative values) . Background events 
are included, and all events are efficiency-corrected. We mea- 
sure Aln£ = -1-19.2 in data (indicated by the vertical line), 
so we accept the spin-1/2 hypothesis. 




-1 -0.8 -0.6 -0.4 -0.2 -0 0.2 0.4 0.6 0.8 1 

cos6^ 



FIG. 10: The helicity angle distribution for 17^(2455)° candi- 
dates in data. The points correspond to efficiency-corrected 
B~ — > X'c(2455)''p candidates. The curves for the spin-1/2 
(solid line) and spin-3/2 (dashed line) hypotheses are over- 
laid. 



CONCLUSION 

We have presented branching fraction measurements 
for the decays B^ -^ A^p and B^ -^ yl+pvr^: 

B(B° -^ Atp) - 

(1.89 ± 0.21 ± 0.06 ± 0.49) x 10~^ 
S(B^ ^ Atp-K^) = 

(3.38 ± 0.12 ± 0.12 ± 0.88) X 10"'', (21) 

where the uncertainties are statistical, systematic, and 
due to the uncertainty in B{A'^ -^ pK'iT^), respectively. 
These measurements are based on 383 million BB events 
produced by the SLAC B Factory and recorded by the 
BaBAR detector. 

If we combine the statistical and systematic uncertain- 
ties only, we obtain B(B° -^ A+p) = (1.9 ± 0.2) x 10-^ 
which is consistent with a previous measurement by the 
Belle Collaboration of B{B° -^ A+p) = (2.2 ± 0.6) x 
10"^ [30]. Both measurements use the same value for 



BiAi 



pK 11+) . However, our measurement for the 



three-body mode, B{B- -^ A+pi:') = (3.4 ±0.2) x lO"'', 
is significantly larger (by about 4cr) than the previous 
measurement from Belle 6(5" -^ A+pTT~) = (2.1±0.3)x 
10""* [25]. The Belle Collaboration measurement uses six 
coarse regions across the B~ -^ A+pn^ Dalitz plane to 
correct for variations in efficiency; we use much finer re- 
gions and see significant variation near the edges of the 
Dalitz plane. This difference in efficiency treatment may 
account for some of the discrepancy between the two re- 
sults. 

One of the main motivations for studying baryonic 
_B-meson decays is to gain knowledge about baryon- 
antibaryon production in meson decays. We have mea- 
sured the ratio of the two branching fractions, 



B{B- -^ A+PTT-) 
B{BO ^ A+p) 



= 15.4±1.8±0.3. 



(22) 



In this quantity the 26% uncertainty in B{A+ -^ pK^Ti+) 
cancels in the branching ratio. 

We have also measured the fractions of B" -^ A+pTT~ 
decays that proceed through a Sc resonance: 



B{B~ ^ i;c(2455)°p) 

B{B- ^A+jm-) 
B{B- -^ 1:^(2800)"^) 



B{B- 



Atpn- 



(12.3 ±1.2 ±0.8) X 10" 



(11.7±2.3±2.4) X 10" 



(23) 



Assuming no interference with direct decay to A+pTT~ , 
about 1/4 of B~ -^ A+pTT~ decays proceed through a Sc 
resonance. 

The order of magnitude difference between the decay 
rates of B^ — > A+pn^ and two-body decays such as 
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5° -^ A+p, B- -^ i:c(2455)°p, and B" -^ ^^(2800)°^ 
is consistent with the theoretical description [16] that 
baryonic B decays are favored when the baryon and an- 
tibaryon are close together in phase space. This inter- 
pretation is also supported by the observation of the 
enhancement in rate when mA^p is near threshold. Al- 
though the A^p threshold enhancement alone could in- 
dicate a resonance below threshold, enhancements have 
been observed in other baryon-antibaryon systems and 
in decays such as e^e~ -^ pjry [1, 2]. Therefore the body 
of measurements indicates that we are observing a phe- 
nomenon that is common to baryon production from me- 
son decays, and possibly common to baryon production 
in general. 

We have used the angular distribution of the decay 
B~ -^ Z'c(2455)°p to study the spin of the S^ baryon. 
The helicity angle distribution is consistent with being 
uniform, which indicates that the S^ has J = 1/2 as- 
suming the ground state /!+ also has J — 1/2 and ex- 
cludes the J ~ 3/2 hypothesis at the > 4cr level. This is 
consistent with quark model expectations for the lowest 
Ec baryon state. 

We also observe an excited Sc state in i?-meson de- 
cays to A'^pTT~ . We measure the mass of this reso- 
nance to be (2846 ± 8 ± 10) MeV/c^ and the width to 



be 



33 



) MeV. It is possible that this observation 



is a confirmation of a triplet of Z'c(2800) states seen in 
ylji"7r~ continuum production [29]. However, the neutral 
1:^(2800)" has a measured mass of (2802 t j) MeV/c^ and 
width of (61 t II) MeV. The widths of the Z'c(2800) and 
the state observed in B decays are consistent, but the 
masses are 3a apart. If these are indeed the same state, 
then the discrepancy in measured masses needs to be re- 
solved. 

Another possible interpretation is that the excited U^ 
resonance seen in this analysis is not the Sd'^SOO)'^ that 
was previously observed. A clear signal is evident for 
B~ -^ i7c(2455)°p decays, but we do not see any evi- 
dence for the decay B^ -^ Z'c(2520)''p. The absence of 
the decay B^ -^ Z'c(2520)"p is in contrast to a claimed 
2.9(7 signal from an analysis by the Belle Collaboration 
based on 152 million BB events [25]. Also, an examina- 
tion of the B~ -^ A'^pTT~ Dalitz plot shows no evidence 
for the decay B" -^ yl+A(1232)~. The 1:^(2520)° is 
a well-established state, and so is the A(1232)++. Both 
are expected to have J ~ 3/2. The Belle Collaboration 
tentatively identified the Uc{2800)^ as J — 3/2 based on 
the measured mass of the state, while there is weak ev- 
idence that the excited S^ we observe is J = 1/2. It is 
therefore possible that B decays to higher-spin baryons 
are suppressed, perhaps due to the same baryon produc- 
tion mechanisms that suppress two-body baryonic de- 
cays, and that the excited S^ state that we have observed 
is a newly-observed spin-1/2 state. 
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